Wetting of Argon on CO2 
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We have studied the wetting transition of Ar adsorbed on solid CO2 by means of high- precision 
adsorption isotherms measured with a quartz microbalance. We observe triple-point wetting. At 
variance with many theoretical studies based on a model adsorption potential, which predict for 
this system a genuine prewetting transition around 100A", we find that a detailed density-functional 
calculation employing a more realistic adsorption potential leads to triple-point wetting of Ar on 
CO2, in good agreement with the experiment. 



PACS numbers: 68.45.Gd;68.35.Rh 

Recently, there has been a tremendous progress in the 
experimental investigations of wetting phenomena. For 
the first time, conclusive evidence for the first-order na- 
ture of the wetting transition, as originally predicted by 
Cahn |iJ long time ago, together with its accompany- 
ing prewetting jumps away from the coexistence line, has 
been obtained in several different systems: quantum liq- 
uid films physisorbed on heavy alkali metals Q , complex 
organic liquids near critical liquid mercury || and 
binary liquid crystal mixtures 

The first evidence of such a transition came however 
from mean-field numerical calculations of Ar adsorbed 
on solid CO2 carried out by Ebner and Saam || . Several 
groups have subsequently studied this same model sys- 
tem with different numerical and analytical techniques 
and have essentially confirmed the original picture . 
The most recent and accurate calculation predicts a 
prewetting transition of Ar on CO2 around 105° K ||. 
It is curious that despite such a keen theoretical inter- 
est, the wetting properties of Ar on solid CO2 have never 
been investigated in a real experiment, to the best of our 
knowledge. 

We have thus decided to study such a system with a 
quartz microbalance technique. This is a powerful quan- 
titative probe in which the resonance frequency of an 
AT-cut crystal excited in its shear mode depends, among 
other things, on the mass of the film adsorbed on the 
metal electrodes of the quartz plate. To drive the crystal 
to its resonance frequency we have used an FM-techniquc 
described in detail elsewhere H. The quartz plate em- 
ployed in this study has a fundamental resonance fre- 
quency of 5MHz with optically polished gold electrodes. 
The frequency shift caused by the adsorption of a mono- 
layer of liquid Ar, — A/ m , corresponds to approximately 
6Hz. The resolution is of about ±0.05Hz with an excita- 
tion power dissipated onto the crystal of less than 50nW. 
The AT-plate is housed in a double-wall OFHC copper 
cell (see Fig. ([[])) to reduce thermal gradients. The en- 
tire set-up is attached to the cold flange of an homemade 
liquid nitrogen cryostat. The maximum temperature dif- 
ference between the top and bottom parts inside the inner 
copper cell is estimated to be less than 1 [iK. The tem- 
perature stability is better than 0.5 mK. The Ar gas is 
admitted into the cell through a thin capillary very well 



heat sunk to the base of the copper cells. A copper foil, 
glued to the bottom of the inner cell, faces the extremity 
of the Ar capillary to shield the quartz crystal from the 
incoming warm vapor which is slowly admitted into the 
cell during the measurements. 
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FIG. 1. Total film thickness of Ar on CO2 as a function 
of T. The vertical line indicates T t . Inset: schematic drawing 
of the sample cell housing the quartz crystal and enlargement 
of the CO2 evaporator. 

The CO2 substrate is evaporated at low temperature 
(« 85°K ) directly onto the two opposite gold electrodes 
of the quartz crystal through a different capillary also 
connected to the gas system. At the temperature of 
evaporation, the vapor pressure of CO2 is ~ 1/j/Torr. 
CO2 is thus admitted to the cell through a copper cap- 
illary mounted inside a stainless steel tube, soldered to 
the base of the inner cell (see Fig. (]]])). The copper 
capillary, wrapped by a resistive wire, is soldered to the 
opposite end of the s.s. tube . In this way, it is possi- 
ble to heat up the CO2 capillary to a temperature much 
higher than that of the quartz and thus avoid its blocking 
during the evaporation. The capillary ends in two curved 
pieces placed in front of the two faces of the quartz crys- 
tal (see Fig. ([!])). By slowly admitting high-purity CO2 
gas through this heated capillary it is then possible to 
evaporate CO2 films with a mass corresponding to thick- 
nesses ranging from 100 to 1000 layers over a period of 
about 2-10 minutes. 

The wetting transition of Ar on CO2 has been experi- 
mentally investigated by means of adsorption isotherms. 



1 



These are determined by slowly admitting small amounts 
of high-purity Ar gas into the sample cell, kept at a con- 
stant temperature. The resonance frequency and the 
pressure inside the sample cell are continuously moni- 
tored with a personal computer. Saturation is reached 
when P remains unchanged after further admission of a 
small amount of gas. 

Fig.(|l|) shows the results of a set of adsorption 
isotherms taken across the bulk triple-point of Ar on gold 
preplated with more than 100 equivalent layers of CO2. 
The vertical axis reports the total frequency shift — A/ tot 
measured at pressures between vacuum and coexistence 
vapor pressure and corrected by vapor contributions in 
the standard way ||||]. Since the film thickness is pro- 
portional to the corrected frequency shift, the maximum 
thickness of the Ar liquid films is found to remain prac- 
tically constant over the investigated temperature range, 
as the horizontal line clearly shows. Finally, on bare gold 
— A/ tot is only 5% smaller than on this CO2 plating. This 
difference is likely to be caused by the rough CO2 sub- 
strate evaporated at low temperature which makes the 
active surface area of the microbalance larger than that 
of the gold electrodes. However, such a small effect indi- 
cates that the roughness of this CO2 substrate does not 
play a significant role in our measurements. 
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T = 84.30 K 
x = -0.343±0.012 
q= 1.03+0.03 
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FIG. 2. Determination of the growth exponent x of Ar on 
CO2 

In the solid phase, —Aftot is found to grow continu- 
ously as T approaches T t . We have also taken isotherms 
on gold preplated with different CO2 coverages ranging 
from 200 to about 1000 equivalent layers. These mea- 
surements agree well with the results of Fig. (|l|) apart 
from an increased frequency noise and a larger value 
of —Af tot observed at the higher coverages, probably 
caused by rougher surfaces. Similar behavior, seen in all 
physisorbed systems characterized by strong substrates, 
is typical of triple-point wetting. We have thus fitted the 
points below T t in Fig. (Q) according to the power law 
expected for this phenomenon p0|, —Aftot = [AT(T t — 



T) x , where A 



OtAr-C0 2 



dT 



liq 



aAr-co 2 being the adsorbate-substrate Hamacker con- 
stant, —Af m = 6Hz the frequency shift corresponding 
to the adsorption of an Ar monolayer and x — —1/3 for 
van der Waals interactions. The results of the nonlinear 
least square fit are: T t — 83.80 ± 0.02i"T, in very good 
agreement with the tabulated bulk triple-point of Ar, 
T t = 83.806 K , x = -0.337 ± 0.008 and A = 1.5 X 10" 5 
K~ 2 Hz~ 3 . We have also analyzed the film growth of 
the isotherms taken above T t . Under the assumption 
of complete wetting at coexistence, the film thickness 
is predicted to grow according to the Frenkel-Halsey- 

-1/3 



-A/ = -Af„ 



Hill relation ]!(_,. _, _,,„ _ 

where the exponent —1/3 reflects again the van der 
Waals nature of the interaction. In Fig.(||) the quantity 
log 10 (—A/) is plotted as a function of log 10 [\n(P v /P)]. 
As expected, near saturation the experimental points lie 
on a reasonable straight line. The slope of the linear 
fit is -0.343 ± 0.012, in good agreement with the FHH 
equation. From the intercept q we have also estimated 
a value for aAr-co 2 — 19000KA 3 ± 20%, in reasonable 
agreement with the van der Waals tail used by Ebner and 
Saam ||. A similar value, although affected by a larger 
error, has been deduced from the parameter A of the the 
previous fit. 

Motivated by these experimental findings, which are 
in disagreement with the previous theoretical predictions 
of a prewetting transition in this system, we have re- 
examined the basis for such predictions and found that 
they are flawed in one of their basic premises, i.e. in the 
form of the fluid-substrate interaction potential. Com- 
mon to these theoretical studies is in fact the assumption, 
as made by Ebner and Saam || , that the Ar-C02 interac- 
tion is represented by a (12-6) Lennard- Jones potential 
whose hard-core diameter and well depth are deduced, 
through simple combining rules pT| , from those of Ar- 
Ar and CO2-CO2. An integration of the (12-6) potential 
over a continuum of CO2 substrate atoms produces a 
(9-3) adsorption potential, hereafter called Ebner-Saam 
(ES) potential, which has been used in almost all the 
theoretical calculations published so far . 

By comparing this potential with a more realistic one, 
which we compute by explicitly taking into account the 
structure of the CO2 substrate, as shown below, we find 
that the ES potential underestimates the real interaction 
between Ar atoms and solid CO2. Our density func- 
tional calculations encompassing the newly calculated 
adsorption potential support the experimental evidence 
reported above of triple-point wetting. 

In our density functional calculation the free energy of 
the fluid is written in term of the density p(r) of the fluid 
as: 

F[p]=F HS [p] + ( P {?)V s (r)dr 



+ -J J p(f)p(r')u a (\r-r'\)drdr' (1) 
Here Fjjs 1S a non-local free-energy functional for the 



2 



inhomogeneous hard-sphere reference system, V s (r) is the 
external adsorption potential due to the surface, while 
the third term is the usual mean-field approximation for 
the attractive part of the fluid-fluid intermolecular poten- 
tial, u a . Fhs is evaluated according to the theory of Ref. 
IjJI , which has been found to give an accurate description 
of inhomogeneous fluids even in situations where strong 
density variations occur, as in the presence of hard walls 
or strongly attractive substrates. Following Ref. ||[l3| we 
describe the attractive part of the interparticle potential 
as an effective interaction 

u a (r) = , r < A 1/6 cr 

= 4e{A(a/r) 12 - {a/rf} , r > X 1/e a (2) 

We determine the free parameters A and a (the HS 
diameter) for each temperature by requiring that the ex- 
perimental values of liquid and vapor densities at coex- 
istence, pi and p v , are reproduced for the bulk fluid. In 
order to correctly describe the bulk phase diagram, it 
is important to verify that the thermodinamic equilib- 
rium conditions are satisfied. This is done by means of a 
Maxwell (equal-area) construction in the P — p plane. 

CO2 has a unique crystalline phase at low temperature 
and pressure (m). Thin solid films of CO2 are known 
to have crystalline structure jl5f when grown at T > 
80° K and in a wide range of growth rates, including 
those realized in our experiment. Moreover, it is found 
experimentally Jl5| that the exposed surface is almost 
invariably the (001) face of the crystal. 

Accurate ab initio potentials are available [Q, which 
describe the interaction between a CO2 molecule and an 
Ar atom. Due to the nature of the bonding in CO2 
solid (van der Waals + electrostatic forces) no apprecia- 
ble electron charge redistribution is expected between the 
molecules in the solid. We may thus calculate the total 
Ar-surface potential by directly summing the two-body 
interactions between an Ar atom and the CO2 molecules 
in the crystal, assuming that the exposed surface is the 
ideal (001). The resulting potential energy surface as a 
function of the lateral position of the Ar atom exhibits 
large variations, with differences in binding energies for 
different adsorption sites of more than 500 °K. In par- 
ticular, in the c(2 x 2) surface unit cell of the clean (001) 
surface there are two (non-equivalent) adsorption sites 
with very large (~ 1000°^) binding energies. We then 
make the assumption that at the low temperatures in- 
vestigated here, at least one Ar monolayer is adsorbed 
on these sites as a solid-like layer. We determine the Ar 
equilibrium positions by energy optimization, neglecting 
however any relaxation of the CO2 surface. The structure 
of the adsorbed monolayer is shown in the upper part of 
Fig. (||). We finally average laterally the interaction po- 
tential between Ar and this Ar-preplated surface, to get 
the adsorption potential V s (z) to be used in our calcula- 
tions (z is the coordinate normal to the surface plane). 
This is shown in the lower part of Fig. (||), where it is 
compared with the original ES potential. The origin of 



the z axis is taken at the position of the outmost layer 
of C atoms in the case of the ES potential, while it is 
taken at the position of the solid Ar monolayer for our 
calculated V s (z). Note that the net adatom-surface po- 
tential is deeper than the ES potential even if the CO2 
surface is screened by the solid Ar layer. This is a conse- 
quence of the large corrugation of the CO2 surface: the 
structure and orientation of the CO2 molecules on the 
surface plane (see Fig.(^)) allow the Ar monolayer to be 
easily accomodated but still leaving space for additional 
Ar atoms to come close to the surface, thus experiencing 
the attractive part of the interaction. 
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FIG. 3. Upper part: top view of the CO2 surface. The 
left part shows the clean, ideal surface: C atoms and O atoms 
are represented by the green and red balls respectively. The 
right part shows the surface covered with a monolayer of Ar 
atoms (blue balls), in the positions as explained in the text. 
Lower part: fluid-substrate laterally-averaged potential. The 
dashed line shows the ES potential, while the solid line is our 
calculated potential. See text for the definition of the z — 
planes. 

The equilibrium Ar density profile p{z) is determined 
by direct minimization of the functional (nh with respect 
to density variations. In practice we fix the coverage 
r = J dzp(z) and solve iteratively the Euler equation 
fx = SF/Sp(z), where the value of the chemical potential 
p is fixed by T. 

We have considered two isotherms, one at T = 105 °K, 
which has been studied in most detail in Refs. ||], and 
one at T = 85 °K, i.e. just above the triple point T t . We 
report our results in Fig. (Q) , where the chemical poten- 
tial, measured with respect to its value at coexistence, is 
plotted vs. the film thickness. The latter is expressed in 
nominal layers I = p l 2 ^ 3 J^°[p(z) — p v ]dz. The dots show 
the calculated points at T = 105° K, where one observes 
a smooth increase of coverage with chemical potential as 
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the saturated liquid is approached from below. At vari- 
ance with the results of Ref. ||, we do not observe any 
sign of a prewetting transition. We show for compari- 
son in the same figure the results obtained, at the same 
temperature, by using the ES potential (triangles). The 
large oscillation with positive values of fj, is representative 
of a thin film-thick film prewetting transition, in agree- 
ment with previous calculations using the same potential 
||-[|. A jump in the coverage from I ~ 0.8 to I ~ 6 can 
be readily obtained from our data by means of a Maxwell 
equal-area construction (dashed line). 



and Saam by summing isotropic Ar-C02 LJ interactions 
over a continuum substrate. The increased binding of the 
surface changes dramatically the Ar growth mode, which 
now results in a continuous (layer-by-layer) film growth 
rather than exhibiting a prewetting transition. The re- 
sults of our calculations are consistent with a wetting 
temperature close to the triple point, in agreement with 
our quartz microbalance measurements. 

We thank Stefano Sitran and Gilberto Schiavon for the 
preparation of the quartz crystals used in this work. 
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FIG. 4. Calculated adsorption isotherms of Ar in the fluid 
phase. Triangles: T = 105° K, as calculated with the ES po- 
tential. Dots: T = 105°K, with the new potential. Squares: 
T = 85° K, with the new potential. 

The calculated points at T = 85°K are shown with 
squares. A sequence of oscillations occurs in this case, 
whose period is one nominal layer, corresponding to 
the condensation of successive Ar monolayers ("layering" 
transitions). 

Additional theoretical support to triple-point wetting 
of Ar on CO2 can be obtained from a simple but surpris- 
ingly accurate heuristic model jl7j where the energy cost 
of forming a thick film is compared with the benefit due 
to the gas-surface attractive interaction, 



(pi - Pv) / V s (z)dz = -2-i 



(3) 



Here p\ and p v are the densities of the adsorbate liquid 
and vapor at coexistence, 7 is the surface tension of the 
liquid and z$ is the equilibrium distance of the potential 
V s . By using for V s our calculated potential and the 
experimental values for pi(T), p v (T) and j(T), we find 
that Eq.(§) is satisfied at T w = 1.03 T t . 

In conclusions, our experimental results suggest triple- 
point wetting of Ar on solid CO2, in contrast to the many 
previous theoretical investigations. To clarify this find- 
ing, we have implemented density functional calculations. 
A crucial ingredient is a realistic form for the fluid-surface 
potential, which is obtained by explicitly considering the 
microscopic structure of the surface of solid CO2 and 
the anisotropy of the Ar-C02 interaction. The result is 
a stronger potential than originally estimated by Ebner 
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